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ABSTRACT 


A complex atmospheric turbulence model developed by the RAE for use in 
aircraft simulation is analyzed in terms of its temporal, spectral, and 
statistical characteristics. First, a direct comparison is made between 
cases of the RAE model and the more conventional Dryden turbulence model. 
Next the control parameters of the RAE model are systematically varied and 
the effects noted. The RAE model is found to possess a high degree of 
flexibility in its characteristics, but the individual control parameters 
are cross-coupled in terms of their effect on various measures of intensity, 
bandwidth, and probability distribution. While a guide is given to the 
setting of model control parameters, the resulting effect on pilot opinion 
is not known. 
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SECTION I 


INTRODUCTION 


A complex atmospheric turbulence model has been developed by the Royal 
Aeronautical Establishment (RAE) for use in manned and unmanned aircraft 
simulation. This model, which is described in Ref. 1, has been used to some 
extent in this country and in Europe (e.g., Refs. 2 and 3), Unfortunately, 
using the RAE model involves setting several parameters, the effects of which 
are not well understood. 

The purpose of the study reported here, therefore, was to improve our 
understanding of the characteristics of atmospheric turbulence represented 
by the RAE model and its parameters . This was accomplished by considering 
the model in terms of its temporal, spectral, and statistical characteris- 
tics. A variety of metrics was applied. In addition, a direct comparison 
was made between one particular RAE model formulation and the widely used 
Dryden turbulence model (Ref. 4). 

The comparison between the RAE and Dryden turbulence models is presented 
in Section II with the RAE model parameters set to yield Gaussian-lilce 
characteristics of the turbulence. In addition to considering the two models, 
per se, this section also contains a discussion of the spectral and statis- 
tical parameters used to quantify the characteristics of the turbulence models 

Section III describes in great detail the effects of perturbing each of 
the RAE model input parameters in terms of temporal, spectral, and statis- 
tical parameters used to quantify the characteristics of the turbulence models 
The perturbations of each of the RAE model input parameters are made with 
respect to the baseline settings up'd in Section II. The reader who is pri- 
marily interested in an overview of the effects of the RAE model inputs may 
omit Section III and proceed directly to Section IV. 

Section IV gives a summary of the analyses reported and recommendations 
of how the RAE model should be used in a real-time simulation environment . 
Recommendations for future analyses of atmospheric turbulence models are 
also presented. 
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SECTION II 


COMPARISON OP RAE AND DRYDEN MODEL BASELINE CASES 


Baseline cases fox’ the RAE and Dryden models are compared in this section. 
The main objective is to call attention to major features of the RAE model 
and to present the various means of quantifying those features . The Dryden 
model is used as a point of comparison because it is a familiar standard 
and is relatively simple to quantify. 

The numei'ical values of model parameters were chosen to make the temporal, 
spectral, and statistical characteristics of both model baseline cases ap- 
proximately equivalent . Since the RAE model is fundamentally different from 
the Dryden model in its mathematical form and involves more parameters, true 
equivalence among all the above mentioned charactei'istics is difficult to 
achieve. The following specific conditions were chosen: 

Airspeed, 250 ft/sec 

Altitude, 800 ft 

u-Gust intensity, 1 ft/sec rms 

These numerical values were then used to pick the model parameters given in 
Tables 1 and 2. 

The baseline cases are analyzed in ternis of their temporal, spectral, 
and statistical chai-acteristics . The results are presented in the following 
subsections . 

A. TEMPORAL DESCRIPTION 

Two sets of time histox'ies are given for each of the two baseline cases . 
First, a complete time history of u- and w-gusts is given for the RAE and 
Dryden models . This is followed by a plot which illustrates the effect of 
one of the gust components on a simplified pilot -aircraft model. 
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TAEIE 1 


RAE MODEL BASELINE PARAMETERS* 

Bandwidth control, NGUSTS t = 5 gusts/sec 
Probability distribution control*, F = 0.0 
Decay ratio control, R = 0.7 
Altitude, HCG = 800 ft 
Frame time, DT = 0.1 sec 
Intensity control, SIGMA = 1 .0 


* The model parameters given in this table are the subject of Section III 
and are explained fully therein. 

t A rule of thumb for setting NGTJSTS is: NGUSTS = 

t The term "global intermittency" is applied to this parameter in Ref. 1, 
but, based on our analysis, we consider "probability distribution control" 
to be more descriptive. 


TABLE 2 

DRYDEN MODEL BASELINE PARAMETERS 

5/ g 

Scale lengths, L u = y 1750-800 

■ and L = 800 ft 
w 

Intensities, a = 1 ft/sec and a 

li W 

Spectral forms; 


= 1348.1 ft 


-VF 


= 0.770 ft/sec 



5 
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Cowrote time histories of the two baseline oases two shown in Pi g. 1 . 
'there is no clear Oiatimti«a to to made between the baseline BM3 and Bryden 
models in this figure, Magnitudes and apparent frequency content of u« and 
ww gusts two comparable within each of tho models as well as one model against 
the other, SL'hat is to soy, it would ho difficult to distinguish specific 
gust components ok distinguish tho models themselves from the time histories . 


(Naturally this applies only to the haae3i.no cases chosen here — it is not 
a general observation, ) She spectral and statistical descriptions to follow* 
will batter diatUl th. distinguishing features which do indeed exist. 


In order to obtain a first order sppmimation of the effect of vvgusts 
on an aircraft.,, the following procedure was applied, I'he eloaod loop pilot- 
aircrafti combination was represented, hy a plteh-attitudo-eonatraiwod transfer 
function of the altitude response to a u~gunt, i,e., 



Where " is the partial derivat ive of r - force due to a 
v u-gvwt (approximately -PgA r ) 

sr* ia the low* frequency attitude numerator aero 
■d| (approximately equal to speed damping) 

i 

'XT' ia the high frequency attitude numerator aero 
'dig (approximately equal to heave damping) 


Por the cases presented here* 

If, sj -o «?/aec 

•* 

rs 0.1 rad/ see 

'rt— st 0 . rad/sec 
''•Of; 

Niger© 0 shows a time history of u-gnst and the corresponding altitude 
excursion based on tho above method, Any Gifforeaees between the models 
are not discernible from these time histories. 
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Figure 2. Comparative Time Histories for u-Gust and Altitude Excursion Due to u-Gust 

for RAE and Dryden Model Baseline Cases 



B. BPEOIRAL PEOORIPTION 


The turbulence resulting from the two models were spectrally analysed 
using a fast Fourier transform routine and the results presented in a variety 
of ways. From past experience (Ref. 5 ) we had found that plotting linear 
00 versus log <n provided tetter insights to the spectral distribution 
than did the usual plot of log o(u>) versus log to. This was because the 
former way of plotting showed the true proportion of spectral power versus 
frequency while, at the same time, utilising a log scale for frequency. 
Finally, a third way of presenting spectral information is presented here: 
a plot of cumulative power fraction (i.e., the fraction of the total power) 
Versus frequency as defined below . Cumulative power fraction, PR: 



This is perhaps the most precise way of describing spectral power distribution. 

Figure 3 shows comparative power spectral density plots of log 0 u (tc)/o^ 
versus log <0 for the baseline ME and Bryden cases . Both cases exhibit the 
features of a typical first order low-pass filter, as shown by the solid 
line representing ideal Bryden model power spectra. From this figure it 
is difficult to observe any significant difference in amplitude or break 
frequency between the two b as aline cases . Some disparity appeal's in the 
randomness of the estimated power spectra for u-gust . 

The same observations apply to the w-gust and the coherence between u~ 
and w -gusts (Fig. 3). There is no apparent difference between the two 
models — only the normal amount of randomness in the low frequency spectral 
measurements . 

A better distinction in spectral characteristics is shown by the plots 
of wo(a) versus frequency shown in Fig. 4. Although the dispersion in 
spectral measurements is prominent, it is possible to see more spectral power 
in the vicinity of 1 rad/seo in the MB u-gust baseline case. For w -gusts. 
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Figure 3. Power Spectral Density and Coherence of u- caxd w-Gusts 
Resulting From the PASS end Dryden lorhuleiico Models 
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the same difference in power spectral density 'between the two models appears 
with one exception at 1 .69 rad/sec. 

Our host resolution of measured spectral power distribution came about 
from consideration of cumulative power fraction versus frequency (sea pre- 
vious definition) . The three fractions of the total power to be considered 
were 0.10, 0..5Q, and O.90*. The resulting power fraction plots are presented 
in Fig . 5 for u- and w- gusts . idealised Dryden model characteristics are 
shown by a solid line. 

The spectral power in the BASS model clearly is centered about a higher 
frequency for both gust components . The implication is that a rule of thumb 
for setting the bandwidth control, NGUSTS, i.e., NGUSTS » does not 
provide a half-power point equivalent to the Dryden model under the same 
conditions. In fact, the rule of thumb yields a half -power frequency for 
u-gust which was higher than that for the Dryden model by a factor of two 
in this case. Most important, this is not clearly visible in the traditional 
plot of log t versus log to (Fig. 3) • 

0. STATISTICAL DESCRIPTION 

The most interesting comparisons between the RAE and Dryden models wore 
found in the various statistical descriptions which were considered. 

1. Probability Distribution Function 

The basic u- and w-gust probability distributions for both models are 
shown in Fig. $. The scale of the probability axis was chosen such that a 
normal distribution would plot as a straight line. It is seen from Fig. 6 
that the Dryden model used here is not perfectly Gaussian, although it should 
be. The RAE model deviates further from a Gaussian distribution, but not 
dramatically. We shall see a greater deviation from a Gaussian probability 
distribution for other RAE model parameters considered in Section 111. 


* Due to the coarseness of the spectral estimates, it was not always possible 
to obtain the frequency exactly at those values of the power fractions (e.g., 
the frequency at a power fraction of 0.15 instead of 0.10 would be obtained) . 
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Power Traction of w Variance Power Fraction of u Variance 



D <1 







2 « Time-Stationarity 


One method, used to examine statistical characteristics was to apply a 
"sliding temporal window" to the turbulence samples considered. Thus over 
a window having a duration of T seconds, various statistical relationships 
were measured and expressed as functions of T. 

The sliding window was considered to be one device for measuring the 
essential features connected with wind shear. Windows much shorter or much 
longer than predominant pilot -vehicle response time constants should be of 
lesser consequence than windows which are approximately equal to those time 
constants. Naturally, the precise range of important pilot-vehicle response 
time constants (or frequencies) depends upon the specific aircraft, its 
configuration and flight condition, and, most importantly, upon the pilot* s 
closed loop structure and degree of tightness. If normal landing approach 
flight path and airspeed regulation is involved, the respective crossover 
frequencies will be about 0.1 rad/sec and 0.5 rad/sec. The corresponding 
significant window length, therefore, would presumably be about 5 to 10 sec 
in duration. 

Another way of interpreting window duration is its relationship to the 
average duration of an approach. Suppose a landing approach normally takes 
one to two minutes. We could therefore concern ourselves with statistical 
properties over that period. The above concepts are illustrated in Fig. 7 - 

We shall see that such properties can vary appreciably from "steady- 
state" values . Furthermore, measures of stationarity are provided by the 
mean and dispersion of various statistical measurements (e.g., the variance, 
kurtosis, etc.) due to sampling the turbulence via different window durations. 
Convergence of the mean of the statistical measurements to a limit is one 
indication of stationarity. In addition, a low dispersion of the statistical 
measurements is another indication of stationarity. Reference 6 contains a 
more detailed discussion of conditions for stationarity of stochastic 
processes . 

a. Standard Deviation . Standard deviation was the most elementary 
statistical property measured as a function of window duration. Figure 8 
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Figure 7* Generic Example of Sliding Window' Concept 
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Figure 8. RMS u- and w-Gusts for the Dryden and RAE Turbulence Models 
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shows rms u-gust in terms of its own average (o^) and standard deviation ( crcr u ) 
relative to the overall rms (cr u ) for varying window duration, T. Note that 
as T approaches the total sample duration, T (204.8 sec): 




and 




u 


u 


-► 0 


The sample size limitation causes o - ^ — ► 0 as T — >- i’ run (e.g., the variance 

of one sample is zero) . If T were doubled then a 0u would approach zero as 

T — *► 2 T , but would be non-zero at T = T . On the other hand, for 
run run 

T«T run 'the rms is smaller and the variability is larger. For example: 




a u 


°du 

T 


°u 


a u 


RAE 

Dryden 

RAE 

Dryden 

8 sec 

0.7 

0.6 

0.25 

0.20 

64 sec 

1 .0 

0.9 

8 

• 

o 

0.15 


We note that the Dryden rms ratio is somewhat smaller (i.e., a^/ar is 
smaller) than the RAE rms ratio for any given window length. This is con- 
sistent with the difference in half -power points noted earlier. Also, cr 0u 
for the RAE u-gust is consistently greater than that for the Dryden u-gust, 
which indicates that the RAE turbulence is less stationary than the Dryden 
turbulence . 

Figure 8 shows the rms dependence on window length for w-gust, and we 
can observe less difference between the models than for the u-gust . Again, 
this is consistent with the spectral differences noted earlier . 

b. Kurtosis . Kurtosis, a measure of the fourth order central moment, 
was calculated in terms of the sliding window and was more sensitive to some 
model differences than variance. Figure 9 shows comparative kurtosis for 
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□ RAE Baseline 
A Drydon Baseline 
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Figure 9. Kurtosis of u- and -w-Gusts for the 
Dryden and RA.E Turbulence Models 
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u-gust and vi -gust . Kurtosis for a normally distributed random variable 

* 

should be throe, and this is the ease for u- and v-gust as 1' becomes large 
for the Drydcn turbulence , (Kurtosis for the entire Drydon model samples 
were measured as 5 .1 and 8.9 for u- and w -gusts, respectively , ) fhe mean 
hurt os is of the ME model was found to be larger than that of the Dryden 
model for longer duration windows but not for shorter durations, say, less 
than 16 see. In the latter ease, the mean kurtosis was about equal . As a 
result wo might aspect to find little evidence of so-called non-Gaussian 
effects in the short term compared to the very long term or steady state. 


Another feature plotted in Fig. 9 is the variation of kurtosis as a given 
window is moved along the finite duration samples . DJhls feature is shown by 
the vertical bars extending from the mean kurtosis . The variation in n-guat 
kurtosis becomes very large for short windows in the case of the ME model. 

Wq find it difficult to relate this statistical feature to a specific time 
domain feature but believe that it may reflect the largo amplitude rates 
of change in u-gust which occur from time to time in the ME model. It may 
also relate to the feature to be discussed next . 


5 , Effective Wind Shear 


Wind shear was one characteristic of particular interest in this study. 
In previously unpublished work (Kef, 7) we had derived analytical relation- 
ships for the Drydon model to characterise the effective change in gust 
wi\plitude over a given time interval (window) . liwo ways of expressing this 
feature are illustrated in Fig. 10. SPho first way involves a least squares 
fit of u(t) over a given time interval and thus represents an effective rate 
of change of u-gust . five second metric takes the simple increment in u over 
a given time interval, both metrics are averaged over a finite duration run 
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using the sliding window idea, For an infinitely long sample of a Dry&en 
model u-gust the following relationships can be derived: 


1 . For a gij.'j: 

S ),. o 2 r .. ] 

* 3 ** fl S ) - 3 (a + ii) a e" 1 j 

where i| » 1—^1 

rvl 

2. For Abg w n(tHl') - u( t ) 



So u S (1 - o-») 


Figure 11 shows the first and second parameters plotted as a function of 
1’ for the ME and Di'yden turbulence models . Ewe opt at small values of 1’, 
neither plot shows any consistent differences between the two turbulence 
models . 


fills concludes the comparison of the Dryden and ME turbulence models 
and the description of the spectral and statistical measures used to quantify 
the characteristics of the turbulence models , The next section describes 


the effects of perturbing each of the ME model input parameters in terms 


of the spectral and statistical measures already described. 


TR- 11 26-1 


20 



□ HAE Bor.olino 
A Prydon Msolina 


10 -.0 
1' (soo) 


-AD 


□ KAK Bivr.Qlino 
A Diryd.au Bytsolina 


1’ (see) 


Figure 1 1 . RMS Velocity Increments for the 
Dryden and RAE Turbulence Models 


TR-USS-I 


21 



SECTION III 


EFFECTS OF VARYING THE RAE TURBULENCE MODEL PARAMETERS 


The effects of individually perturbing each of the RAE model parameters 
on the temporal, spectral, and statistical properties of the turbulence are 
presented and discussed in the following subsections . The perturbations are 
made with respect to the baseline settings defined in Section II, which were 
intentionally chosen to yield Gaussian -like characteristics of the RAE 
turbulence. Thus by perturbing the RAE model parameters with respect to 
these nominal settings the non-Gaussian characteristics of the RAE turbulence 
will be demonstrated and quantified. The nominal and perturbed values of the 
six RAE model parameters are shown in Table 3 . 


TABLE 3 


BASELINE AND PERTURBED VALUES OF THE RAE 


MODEL PARAMETERS 


Parameter 

Baseline Value 

Perturbed Values 

NGUSTS 

3 gusts/sec 

2, 8 gusts/sec 

F 

0.0 

oJi-, 0.8 

R 

0.7 

0.3, 0.9 

I1CG 

800 ft 

200, 2500 ft 

DT 

0.10 sec 

0 . 05 , 0 . 15 , 0.20 sec 

SIGMA 

1 .0 



* The RAE parameter SIGMA is used only to set the variances of the 
turbulence components. Therefore, perturbations were not analysed. 
(See Subsection IX .F . ) 
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A brief description of how each parameter is utilised 'within the BAB model 
is given at the beginning of each of the following subsections . In addition, 
whenever it is appropriate, each parameter is quantitatively described in 
terms of its major impact on the attributes of the PAIS turbulence . For example, 
the parameter WGUSTS primarily affects the frequency content of. the turbulence 
and we have therefore dubbed it the "bandwidth control." HGUSiS does, however. 


have secondary effects 


on the statistical characteristics of the turbulence. 


Except for the subsection on the effects of altitude, HOG, the spectral 
and statistical properties of the turbulence are described only for the u-gust 
component. 1'his is sufficient because above D^OO ft the statistical, proper- 
ties of all throe components of the BAB turbulence are equivalent . The only 
interesting property of the BAB turbulence we discovered by examining the 
other turbulence components is their unusually high coherence (high with 
respect to the Dryden model in which the turbulence components are intention- 
oily designed to be independent) . l’his high coherence property is presented 
and discussed in Subsection XXX .D. 


When generating the BAB turbulence time histories, which were subsequently 
spectrally and statistically analyzed, the turbulence was always sampled at 
the same point in time and, even though only ltd sec time histories are shown, 
the length of the samples was always S04.8 sec (20 ! lB data points with a frame 
time of 0.10 sec). Slightly different absolute values of the spectral and/or 
statistical characteristics would be obtained if the run length or the sampling 
points were changed, but the trends in the turbulence attributes connected 
with the BAB model parameters reported below should not change . 


A. BANDWIDTH COM’ROXi, MUSTS 

The amount of high frequency power, or "bandwidth," of the BAB turbulence 
is primarily controlled by the ROUSTS parameter. ROUSTS performs many 
functions within the BAB turbulence model* however, its primary function is 
to control the number of new low-, mid-, and high-frequency gusts generated 
each second. These gusts are inputs to special. -purpose filters, and the 
outputs of the filters are then used to form the three turbulence components: 
u-gust, v-gust, and w-gust . ROUSTS can take on only integer values and is 
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limited to 


1 < NGUSTS < l/DT 

where TO is the frame time. The effects of NGUSTS on the temporal, spectral, 
and statistical properties of the RAE turbulence are presented and discussed 
below. 

vi and w-gust time histories for NGUSTS = 2, 5, and 8 gusts/sec are presented 
in Fig. 12. These show that increasing NGUSTS does indeed increase the band- 
width of the turbulence. This observation is confirmed by examining the plots 
of power spectra and power fraction plotted versus frequency. An increase in 
the power at higher frequencies is discernible from the power spectral density 
plots, <Ii uu (a>)/ cT and oa3> /c u versus ui, in Figs. 13 and 14. The effect on spectral 
power distribution is best depicted by the power fraction plot in Fig. 15. 

Note that the frequency at a power fraction of 0.90 of the total variance, 
co(0.$»), monotonically increases from 1.1 to 4.0 rad/sec as NGUSTS is increased 
from 2 to 8 gusts/sec. The shift in the power fraction is not as dramatic 
or consistent at lower frequencies (01(0.15) increased from 0.06 to 0.12 rad/ 
sec for NGUSTS - 5 to 8 gusts/sec, but remained at 0.u6 rad/sec for NGUSTS = 

2 gusts/sec). This indicates that increasing NGUSTS only extends the band- 
width of the turbulence, rather than shifting the entire spectrum to a higher 
frequency (i.e., the low frequency power remains fairly constant as NGUSTS 
is increased) . 

Figure 16 shows that NGUSTS did not affect the probability distribution 
function (PDF) in a consistent manner . That is, the tails of the PDF were 
altered by NGUSTS, but the low probability tail was not changed in the same 
manner as the high probability tail. Also, neither tail was changed mono- 
tonically by NGUSTS . Thus it appears that NGUSTS has an effect on the 
probability distribution functions, and hence will have an effect on the 
higher order moments . 

Figure 17 shows that the mean limited window standard deviation of u-gnst, 
o u , monotonically increases with NGUSTS. (-The change in "o u d,ue to NGUSTS 
could be compensated for by using the intensity control parameter, SIGMA, 
as discussed in Subsection III.F.) The dispersion of the limited window 
standard deviation of u-gust, o Gu , increases slightly as NGUSTS is increased 
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Figure 16. Effects of NGUSTS on the Probability Distribution Function 
of u-Gust of the RAE Turbulence Model 
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Figure 17. Effects of RGU3TS on tlie rms u-Gust of the RAE Turbulence Model 



from 2 to 5 gust s/s cc, but is not noticeably altered as NGUSTS is further 
increased from 5 to 8 gusts/sec. Thus there appears to he an upper bound on 
the level of limited sample dispersion in a u due to NGUSTS . 

The mean and standard deviation of the hurt os is, ct^ and , were not 
consistently or significantly altered by NGUSTS, as is demonstrated in Fig. 18. 

The standard deviation of the effective wind shear, oau^ is presented 
in Fig. 19. It is seen that o^ui was consistently increased as NGUSTS was 
increased from 2 to 9 gusts/sec, but that cr^ did not change consistently 
as NGUSTS was increased from 5 to 8 gusts/soc (cau-j ws increased at some 
values of T and decreased at others). 

B. PROBABILITY DISTRIBUTION FUNCTION CONTROL, F 

The purpose of the BAE modol parameter F is to alter the shape of the 
probability distribution function (PDF) (although it is dubbed the "global 
intormittency parameter" in Ref. 1). As F is increased from its lower limit 
of 0.0 to its upper limit of 1 .0 the shape of the PDF changes from nearly- 
Gaussian to highly non-Gaussian. The effect is to attenuate the magnitudes 
of the high probability regions of the PDF and amplify the magnitudes of the 
low probability regions (i.e., the tails). F also affects, to a certain 
degree, some of the other characteristics of the turbulence, which are dis- 
cussed below. 

Figure 20 contains time histories of u-gust and w-gust for throe different 
values of F: 0.0, 0 Jl , and 0.8. The time histories show that the major 

effect of F is to eliminate many of the low-amplitude gusts and increase the 
large-amplitude ones, which is the intended purpose of the F parameter. 

F appears to have some small and subtle effects on the spectral charac- 
teristics as evidenced by the power spectra and power fraction plots presented 
in Figs. 21, 22, aud 2$ . It is difficult to discern trends from the plots 
of 4> (to) / a or a>T> (to)/ o'" versus to, however, the power fraction plot of 
Fig. 25 shows a consistent trend, albeit small, for the low-frequency power 
to decrease as F is increased . For example, the frequency at a power fraction 
of 0.50 of the total variance, to(0.50), is decreased from O.k to 0.25 rad/ 
sec as l 5 ' is increased from 0.0 to 0.8. The high frequency power, as reflected 
by to(0.90), is unaffected by F. 
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Figure 23. Effects of F on the Power Fraction of the u-Gust Variance 
Versus Frequency of the BAE Turbulence Model 


Figure 24 shows that F affects the PDF as expected; that is, increasing 
F causes the PDF to become less Gaussian-like, as evidenced by comparing the 
PDFs shown in Fig. 2k to a straight line. The probability distribution 
functions show that as F is increased the high probability gusts are reduced 
and the low probability gusts are increased. This is consistent with the 
observations made of the u-gust time histories shown in Fig. 20. Because 
of these non-Gaussian characteristics the higher order moments (e.g., the 
kurtosis) should also be non-Gaussian. 

Figure 25 shows that the limited window standard deviation of u-gust, 

? u , is virtually unaffected by F, but that the dispersion of the limited 
window standard deviation of u-gust, cr CT , is monotonically increased as F 
is increased from 0.0 to 0.8. Furthermore, note that o for any one value 

of F is unchanged for values of T between k and 6k sea. 

Figure 26 shows that the limited window mean kurtosis of u-gust, ch , and 
the standard deviation of the kurtosis of u-gust, cf^ , are only slightly 
increased as F is increased from 0.0 to 0.4, but are dramatically increased 
as F is increased from 0.4 to 0.8. Thus F is a highly non-linear control 
with respect to the kurtosis. It is seen that will vary between about 2.4 
and 9.0 for F = 0.8, but only between about 2.2 and 4.4 for F - 0 .0 or 0.4 
for a window length of 52 sec. 

F has virtually no affect on the effective wind shear parameter, oAu 1 > 
as shown in Fig . 27 • 

0. DECAY RATIO CONTROL, R 

The decay ratio, R, is a parameter in the filter of the RAE model, and 
is used recursively to restore the high-, mid-, and low-frequency components 
of u-gust, v-gust, and w-gust to zero (refer to Ref. 1 for a more detailed 
discussion) . As will be demonstrated below, R affects not only the temporal 
and spectral characteristics but also the statistical characteristics . 

Because of its multi-dimensional effects on the turbulence characteristics, 
it is not recommended that R be varied unless the desired spectral and sta- 
tistical characteristics cannot be obtained with a combination of NGUSTS 
and F. R is limited to values greater than 0.0 and less than 1 .0. 
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Figure 24. Effects of F on the Probability Distribution Function 
of u-Gust of the RAE Turbulence Model 
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Figure 28 contains time histories of u-gust and w-gust for three different 
values of R: 0 . 5 , 0.7, and 0 . 9 . We see that R has a significant effect on 
the magnitude of the turbulence; however, it is difficult to ascertain any 
changes in the spectral or statistical properties of the turbulence solely 
from time histories . The effects of R on the spectral characteristics of 
u-gust are presented in Figs. 29, 30 f and 31 • Rote that increasing R has 
the effect of reducing the amount of high-frequency power. This is best 
demonstrated in the power fraction plots of Fig. 31 . While the effect of 
R on the power fraction of the u variance appears to be non-linear with fre- 
quency, this is not the case. If, for example, the frequency at a power 
fraction of 0.90, n>(0.90), is plotted versus R using linear scales, it can 
be seen that 0 ( 0 . 90 ) is directly proportional to R. The point is that fre- 
quency content, as reflected by the power fractions, can be nearly linearly 
controlled via the decay ratio R. However, as mentioned before, R will also 
affect the statistical properties of the turbulence. 

The effects of R on the probability density function (PDF) are depicted 
in Fig. 32. Changing R from 0.7 to 0.5 does not have a very noticeable 
effect on the PDF, but going from 0.7 to 0.9 does have a definite effect on 
the low probability regions of the PDF (i.e., the "tails"). Thus we can 
expect R to have an effect on the higher order moments of the turbulence . 

The effects of R on the limited window standard deviation, a , and the 
dispersion of the standard deviation, cr au , of u-gust are shown in Fig. 35* 

Both cr u and a a ■ monotonically increase with increasing R, with the biggest 
increment in both parameters occurring for R = 0.7 to 0 . 9 . The parameter 
o u can be linearly altered with the intensity control parameter, SIGMA. 

(see Section III.F), which as no effect on a au . 

The effects of R on t.h. ,■ mean limited window kurtosis, <x^, and the standard 
deviation of the limited window kurtosis, 0 ^, of u-gust are shown in Fig. 34 • 
is monotonically decreased with increasing R at low values of T (a^ =2.7 
to 1 .7 for R = 0.5 to 0.9, respectively at T = 4 .0 sec), but is not discern- 
ibly affected by R at high values of T (e .g ., = 3 *5 to 3 *4 at T = 64 sec). 

0 ^ is not noticeably or consistently altered by R. 
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Figure 30 . Effects of R on a>I> uu (ao) of the RAE Turbulence Model 
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Figure 31 . Effects of R on the Power Fraction of the u-Gust Variance 
Versus Frequency of the RAE Turbulence Model 








Model 


R has a large effect on the standard deviation of the effective wind shear, 

, as is demonstrated in Fig. 35- The effect is non-linear because the 
change in at all values of T was much larger for changing R from 0.7 to 
0.9 than it was for changing R from 0.7 to 0-5* 

D. ALTITUDE CONTROL, HOG 

The altitude parameter, HCG, affects only the w-gust component. The 
low-, mid-, and high-frequency components of w-gust are linearly decreased to 
zero as HCG goes from 2500 ft to 800 ft, from 800 ft to 200 ft, and from 
200 ft to 5° ft above ground level, respectively. Thus w-gust is unaffected 
by HCG above 2500 ft, and is zero when HCG is below 50 ft. The temporal 
effects of HCG on w-gust are demonstrated in Fig. 36. The time history of 
u-gust is included in Fig. 36 for comparison with the w-gust time history for 
HCG above 2500 ft. More about the comparison between these two time histories 
will be said later. 

From the time histories of Fig. 36, one can see that the main effect of 
decreasing HCG is to reduce the amount of low-frequency power contained in 
w-gust. This is confirmed by the power spectral density plots of Figs. 37 
and 39 and the power fraction plots of w-gust versus frequency shown in 

Fig. 39- 

Figure 1|0 shows that decreasing altitude tends to increase the normalized 
amplitude of low probability gusts. Thus, w-gust will become less Gaussian 
with decreasing altitude (in contradistinction to Dryden model turbulence 
where w-gust is Gaussian at all altitudes). 

Decreasing HCG also decreases the total amount of power in w-gust, which 
is demonstrated in the plot of the mean limited window standard deviation, 
o , and the dispersion of the limited window standard deviation, 0^, of 
w-gust shown in Fig. 41 . 

Figure k2 shows that the mean limited window kurtosis, 5^, and the 
standard deviation of the limited window kurtosis, CT a ^ , of w-gust generally 
increase with decreasing altitude. This is true only for T much less than the 
total run length. 
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Figure 35* Effects of R on the u-Gust Effective Wind Shear Parameter 

of the RAE Turbulence Model 
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Figure 39 . Effects of HCG on the Power Fraction of the w-Gust Variance 
Versus Frequency of the RA.E Turbulence Model 
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Figure 40 . Effects of HGG on the Probability Distribution Function 
of w-Gust of the RAE Turbulence Model 
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Figure h2. Effects of ECG on the v-Gvst Kurtosis of the 



Turbulence Ifoael 


Figure l | o shows that altitude has a large and consistent effect on the 
effective wind shear parameter of w-gust, o. . Note that o^ w is monotonically 
decreased vrith decreasing altitude. 

A close comparison of the u- and w-gust time histories shown in Fig. % 
reveals a strange feature of the RAF turbulence model, 'fhere appears to be 
very high correlation between u-gust and w-gust, especially for ICG > 2200 ft. 
Notice that many of the large u- and w- gusts occur at the same point in 
time, have the same shape, but occasionally ore of opposite polarity (i.e., 
simply l8o degrees out of phase). We say that this feature is strange because 
such high correlation between the txu'bulence components are very rare in the 
Dryden model where the components are. intentionally designed to be independent. 

In order to quantify the correlation described above the coherence 
2 P 

function between u and w~ was calculated . * 
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The results for HCG = 2J00, 800, and 200 ft are plotted in Fig, Mb The 

plots of P n p v o(co) versus to contain a lot of scatter, but generally show an 

increase in the coherence as altitude is increased. The best measure of the 

o 

u- and w-gust correlation is reflected by the overall coherence, p ^ , wMch 
monotonically increases with, increasing altitude . 

As a basis for comparison the coherence for a high altitude Dryden model 
turbulence run was calculated (i.e., above 1750 ft so that the rms of u- and 
w-gust would be about equal). The results are presented and compared to the 


* The coherence of u and w did not demonstrate the high correlation because 
of the occasional 180 degree phase shift between u and w noted above. 
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Figure 43. Effects of HCG on the w-Gust Effective Wind Shear Parameter 

of the ME Turbulence Model 
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RAE coherence in Fig. 45* Note that for the RAE model is much higher than for 
the Dryden model. 

The reason that so much has been said about the coherence is because it 
is generally overlooked as a turbulence feature (as mentioned before, the 
coherence of the u- and w-gust of the Dryden model is theoretically zero), 
yet could be an important parameter when evaluating a computer- generated 
turbulence model or analyzing real-world turbulence. For example, the 
aircraft excursions resulting from highly coherent turbulence should be 
larger than those resulting from highly incoherent turbulence because the 
aircraft is being disturbed in all three axes at the same time. 

E. FRAMETIME CONTROL, DT 

The frametime, DT, has strange and totally unexpected effects on the 
temporal and spectral characteristics of the RAE turbulence. Unlike the 
Dryden turbulence model (where the frametime will not affect the turbulence 
characteristics until the break frequency of the shaping filters approaches 
the Nyquist frequency), the frametime in the RAE turbulence model will affect 
the characteristics of the gust characteristics if l/DT is not an integer. 

Tlois phenomenon is demonstrated in the time histories and power spectral plots 
shown in Fig. 46. 

Visual inspection of Fig. 4-6 shows that the time histories for DT = 0.10, 
0.20, and 0 .05 seconds are virtually identical. But the time history for 
DT = 0.15 seconds appears to be shifted in time and there is more high- 
frequency content. This effect is difficult to see from the power spectral 
plots shown in Figs. 4-7 and 48, but is confirmed by the power fraction 
versus frequency plot shown in Fig. 49- Note that the frequency at a power 
fraction of O.90 was increased by about 1.0 rad/sec for DT = 0.1 5 sec but was 
virtually unchanged for DT = 0.05 or 0.20 sec as compared with DT ~ 0.10 sec. 

The low-frequency end of the power fraction curve is about the same for all 
the frametimes . The frametime did not have a noticeable effect on any of 
the statistical properties of the txu’bulence. 

The reason for the strange effect noted above is found in the way the RAE 
model is implemented. The "rise time" of the high-frequency gust component 
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Figure 48. Effects of DT on <uj> (to) of the RAE Turbulence Model 
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is expressed as an integer number of samples of the frametime, 


MRT3 = ( 5/NGUSTS )/DT 


Because of truncation, the true rise time, TRT3 = MRT3 • DT, is incorrect 
for non- integer values of l/DT. For example, for NGUSTS = 5 and DT = 0.10, 
0.20, and 0.15 sec T RT j is ; 


TRT3 


1.0 sec for DT = 0.10 sec 
1 .0 sec for DT = 0.20 sec 
0.90 sec for DT = 0.1 5 sec 


IF. INTENSITY CONTROL, SIGMA. 

The overall intensity is controlled by the parameter SIGMA. The 
variances of the three gusts components, u-gust, v-gust, and w-gust, all 
scale linearly with SIGMA, regardless of the values of NGUSTS, F, R, HCG, 
or DT. None of the other statistical, temporal, or spectral properties of 
the turbulence are affected by SIGMA. It should be emphasized that, unlike 
the Dryden model, the standard deviation of the RAE model turbulence is 
not equal but only proportional to SIGMA . 
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SECTION IV 


SUMMARY AND RECOMMENDATIONS 


The characteristics of the complex RAE atmospheric turbulence model have 
been analyzed and presented in the foregoing sections with two main objectives . 
The first was to compare the RAE model with a relatively well known form, 
the Dryden model. The second objective was then to vary systematically the 
RAE model parameters in order to observe their overall effect on temporal, 
spectral, and statistical properties. It is important to recognize, however, 
that although the numerical properties are analyzed in detail, we still do 
not have a clear definition of how they impact pilot opinion. 

The results of comparing the Dryden and RAE models reveal that the 

\ 

spectral and statistical properties of the individual turbulence components 
can be made to appear similar if the input parameters of the RAE model are 
properly chosen. One striking difference between the two models is the high 
coherence among the turbulence components of the RAE model . The high co- 
herence is visible by examining the time histories of u- and w-gust (or u- 

and v-gust) and is best quantified by computing the coherence between u and 
2 

w . This feature of the RAE turbulence is emphasized because the individual 
components of Dryden turbulence are intentionally designed to lack coherence. 
The high coherence among the components of the RAE model turbulence could 
be an important feature when using it in aircraft simulation investigations . 

The qualitative effect of each of the six RAE model parameters on 
spectral and statistical characteristics are summarized in Table 4. Quanti- 
tative effects are given in Section III . 

The four most influential RAE model parameters are SIGMA, NGUSTS, F, 
and R. SIGMA varies intensity directly but affects no other property of the 
turbulence . NGUSTS varies bandwidth directly and F varies the probability 
distribution directly. R is a parameter which has strong but diverse effects 
on a number of properties and should be adjusted with considerable caution. 
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TABLE 4 


SUMMARY OF THE SPECTRAL AND STATISTICAL EFFECTS 
OF THE RAE TURBULENCE MODEL PARAMETERS 


RAE Model Parameter 


Turbulence Property 


Remarks : 


Spectral: 

High frequency Power 
Low Frequency Pow.er 

Statistical: 

Variance 
Kurtosis 
Stationarity(3 ) 


Bandwidth 

Control 

HGUSTS 


Primary determinant 
of bandwidth hut 
also has a 
noticeable effect 
on intensity. 


Increase^ 2 ) 


Ho change 


Increase 
Ho change 
Ho change 


Probability 

Distribution 

Control 


Primary determinant 
of probability 
distribution but 
has an indirect 
effect on spectral 
properties . 


Ho change 
Decrease 


Ho change 

Increase 

Decrease 


Decay 

Ratio 

Control 

R 

Affects many 
properties 


Intensity 

Control 

SIGMA 


Varies only the 
intensity but its 


simultaneously numerical value 


and should be 
varied with 
great caution. 


Decrease 

Decrease 


Increase 

Decrease 


Altitude 

Control 

HCG 


Affects w-gust 
only but impacts 
many of its 


does not necessarily properties . 
correspond to 
actual rms. 


Ho change 
Ho change 


Increase 
Ho change 
Ho change 


Decrease 

Decrease 


Increase 

Decrease 

Increase 


Frametime 

Control 

DT 


Ho effect 
on turbulence 
properties 
unless integer 
condition is 
not met ( 1 ) . 


Conditional 
Ho change 


Ho change 
Ho change 
Ho change 


Si 

o. 

& ^ 
Q. nj 


0) The frametime should he set such that l/flT is an integer, otherwise the high frequency spectral characteristics are affected. 

(2) The word "increase" refers to the effects of increasing the magnitude of the indicated RAE model parameter on spectral and statistical 
properties . numerical values for the amount of change can be obtained from Section III , 

(3) Stationarity Is reflected by the dispersion of the statistics in sampling over a limited window. Small dispersion implies high stationarity. 



Parameters having effects of lesser importance are HCG and DT. HCG, 
the altitude parameter, mainly varies the intensity and frequency content of 
w-gust. DT sets the turbulence frametime and a number of parameters internal 
to the RAE model. 

Practically speaking, we should point out that the values of RAE model 
parameters are, in reality, restricted. For example, HCG is determined by 
the particular aircraft altitude. Frametime, DT, must not only be compatible 
with the simulation computer frametime but must also be set such that l/DT is 
an integer. From our study, the decay ratio, R, affects all turbulence 
properties; therefore, we recommend that it be set at a fixed value (about 
0/f) and not be used as a "control." NGUSTS, the bandwidth control, should 
be set according to airspeed in order to provide a desired scale length . 

(As in the Dryden model, scale length, hence NGUSTS, could also be varied 
as a function of altitude.) 

In view of the above parameter constraints, we are left with two primary 
controls for the RAE model, SIGMA and F. SIGMA, the intensity control, is 
perhaps the most fundamental determinant of turbulence effects . Certainly 
for the Dryden model its influence has been shown to be paramount (Ref, 5 ) . 

The probability distribution control, F, adjusts the non-Gaussian qualities — 
the relative frequency of occurrence of large magnitude gusts . It too has 
been judged important (Ref. 8 ), but its desired numerical value is lass 
well defined. Several research studies have shown that the non-Gaussian 
characteristics of turbulence can be quite important in a piloted simulation 
(Refs. 8 , 9 , 10, and 11). 

With only two input variables to be concerned with, an experiment could 
be carried out to evaluate the RAE turbulence model in a piloted simulator. 

We recommend that such an experiment be designed to compare the RAE model 
with some of the other so-called non-Gaussian turbulence models (Refs. 3, 8 , 

9, 1 0, and 1 1 ) in addition to comparing it with the Dryden model . As an 
alternative to using a real-time piloted simulation, an experiment could be 
performed using a realistic pilot model in a non-real-time environment similar 
to that described in Ref. 12. This would be cheaper than a real-time piloted 
simulation, and it would permit a researcher to isolate better the effects 
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of the turbulence model without having to cope with unintentional changes 
in piloting techniques , 

Ono additional recommendation is to decouple the effects of the ME 
parameters on the turbulence characteristics, especially those of NGUSTS. 
Note from Table t that changing NGUSTS Increases the liigh-frcquency power 
and increases the variance of the turbulence . Using the results of Sec- 
tion XIX, the RAE model could he programmed in such a way that the variance 
would he unaffected by changes in NGUSTS . Also note from Table 4 that F 
affects the low frequency power as well as the leurtosis and stationer! ty of 
the turbulence. However, the magnitude of the low-frequency power changes 
due to changing F arc fairly small (see Section 111), and hence we do not 
recommend that any model changes he made in order to eliminate this effect . 
It is not known if it is possible (or indeed even desirable) to separate 
the simultaneous changes in kurtosis and stationarity, although further 
research of these characteristics of atmospheric turbulence should he per- 
formed. As mentioned before, the intensity control parameter, SIGMA, is 
indeed an independent control in that it affects only the turbulence 
component variances . 
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